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Abstract 

Recently, muon production in electron-proton scattering has been suggested as a possible 
candidate reaction for the identification of lepton-flavor violation due to physics beyond the 
Standard Model. Here we point out that the Standard-Model processes e~p — * i^pVpVf, and 
e~p — ► e~nii r Va can cloud potential beyond-the-Standard-Model signals in ep collisions. We 
find that Standard-Model ep — > fiX cross sections exceed those from lepton-flavor-violating 
operators by several orders of magnitude. We also discuss the possibility of using a nuclear 
target to enhance the ep — ► fj,X signal. 
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I. INTRODUCTION 



A number of experiments over the past decade provide compelling evidence that the 
neutrino mass matrix is non-diagonal in the basis of weak eigenstates \v a ); a — e, fi, and r 
(see |l| for a recent review). This knowledge has led to renewed interest in lepton- flavor 
violation (LFV), which can be probed by searches for rare decays such as /i — > cy. Such 
LFV decays are possible when the Standard Model is extended to include neutrino mass 
and neutrino mixing, but the resulting cross section is exceedingly small (branching ratio, 
BR ~ 10~ 60 ) as the process scales with the fourth power of the ratio of the neutrino mass 

results from the Minimal Super-Symmetric extension of the Standard Model (MSSM) [3|. 
The MEG (mu->e gamma) experiment at the Paul Scherer Institute (PSI) [4j is capable 
of detecting branching ratios as small as 10~ 13 at a 90% confidence level, and will search 
for the LFV decay [i — > ej. Such searches for lepton-flavor violation potentially offer an 
intriguing window on beyond-the-Standard-Model (BSM) physics. 

Recently, the possibility of observing lepton-flavor violation in fixed-target electron 
scattering has been raised as an alternative j^j to experiments searching for the rare 
p, — > e7 decay. Facilities with electron beams of high intensity and significant duty factor, 
such as Jefferson Lab, seem to be natural places to perform experiments to search for 
ep — > up. Hereafter, we will refer to this electron-to- mu on conversion process as EMU. 

LI " — 

However, the conclusion of Ref. [2j is that even under the most favorable dynamical 
scenario (a heavy right-handed Majorana neutrino with m u ~ 0(miv)), the Standard- 
Model supplemented by dynamics that results in neutrino oscillations yields a cross section 
a ~ 10~ 27 femtobarns (fb) for EMU — so low as to be inaccessible to current experiments. 

In this paper, we discuss two Standard- Mo del processes that can cloud an EMU signal 
in e~p scattering by generating final states fi~X other than the desired final state fi~p. 
In Standard-Model mechanisms, the additional particles X must have baryon number 1, 
muon lepton number -1, and electron lepton number +1. Two such reactions are: (i) 
e~p — > fi~v^i> e p, and (ii) e~p — > e~n^i + v^. From now on, we refer to these reactions as 
"muon-pro duct ion processes" in order to distinguish them from EMU. 

The first reaction involves only electroweak interactions, and takes place as the electron 
goes off-shell in the scattering event by an amount corresponding to the momentum of 
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the virtual photon exchanged with the target. The electron then decays via the weak 
interaction to a muon, accompanied by the emission of v e and V^. We note that this 
process can also take place off a neutron, although in practice this implies a nuclear 
target. In fact, for the nuclear-target case coherent electron interactions with the total 
nuclear charge can enhance the signal. 

The second muon-production process involves the strong interaction. Even when the 
electron energy is below the pion-production threshold the exchanged virtual photon can 
interact with the "cloud" of virtual pions that surrounds the nucleon. This can generate 
an off-shell 7r + , which decays to a /i + and v^. (In electron-nucleus scattering, the presence 
of neutrons allows e~n — > e'p^v^ to occur through virtual 7r~'s.) As the energy of 
the incident electron approaches the pion threshold, the time for which the virtual pion 
lives (and hence the distance it travels before decaying) increases. Consequently, these 
muon-production reactions switch over to pion electro-production at electron energies of 
about 140 MeV swamping any possible EMU signals. 

This paper is organized as follows. Section [III presents a pedagogical calculation of 
electron-proton scattering as a probe of BSM physics using generic low-energy effective 
couplings that can cause EMU. The diagram involving photon exchange plays a dominant 
role, but is severely constrained by the experimental bound on the coupling obtained from 
the ji — > e7 process. (For a specific realization in the MSSM see Ref. [3|.) In Sec. IHH 
we derive the matrix element for the process e~p — > pTv^VeP from electroweak theory, 
compute the size of the cross section, and provide a simple explanation for the order-of- 
magnitude of our result. In Sec. IIVI the virtual-pion production and decay contributions 
to the matrix element for ep — > en^u^, and the resulting differential cross section are 
presented. Due to differences in the interaction couplings and phase-space factors, the 
cross section for e~p — > e~nfi + u^ turns out to be several orders of magnitude larger 
than that for e~p — > {i~v^v e p. Therefore any EMU experiment seeking BSM (or even 
electroweak) physics would either have to veto processes in which a scattered electron 
is detected in coincidence with the produced muon, or, more feasibly, detect the charge 
of any muons produced in the electron-proton collision. In Sec. [V] we describe muon- 
production in electron-nucleus scattering, including the relative importance of collective 
nuclear excitations. We present our summary and conclusions in Sec. IVI1 Details of 
phase-space integrations and numerics are provided in the appendices. 
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FIG. 1: Beyond-the-Standard-Model contribution to muon production ( "electron-muon conver- 
sion"). The hatched vertex represents the dimension- five coupling of Eq. (pQ). The solid lines 
denote leptons, while the double line is the proton. The particle momenta are indicated in 
parentheses. 

II. ELECTRON-MUON CONVERSION VIA PHYSICS BEYOND THE STAN- 
DARD MODEL 

In this section, we consider the differential cross section for the reaction ep — > \ip 
induced by operators which change lepton flavor, and hence are low-energy manifestations 
of physics "beyond the Standard Model" (BSM). All such operators are, by definition, 
dimension five or above, and so they produce cross sections suppressed by (at least) one 
power of m^/A, where A is the scale of the physics that results in lepton-flavor violation. 
We will show that there is a dimension-five operator that could, in principle, produce a 
sizeable ep — > fip cross section. However, in practice, bounds from the non-observation 
of the process \i — > ej preclude any observable muon production via this dimension-five 
BSM operator. 

The BSM operator associated with the decay fi — > c-f can be written as 



where ipi is the lepton field of family I (e or /i), and F a/3 is the electromagnetic field- 
strength tensor. The object v is the Higgs vacuum expectation value, and A is the scale 
of the BSM physics that induces this operator. The Higgs vacuum expectation value 
appears because while the operator is of dimension five, it is suppressed by an additional 
power of v/A because it changes lepton chirality. Dimension-six BSM structures which 
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have the low-energy form 

Cf*^ ~ JLfco^NON , (2) 

where N denotes the nucleon field and the O's are operators (potentially with Lorentz 
indices that are contracted with one another) can also appear in Lj. However, their effects 
are suppressed relative to the operator in Eq. ([I]). 

The Feynman diagram for ep — > fip for the coupling in Eq. (CQ) is shown in Fig. [U The 
general form of the nucleon current (given parity invariance, time-reversal invariance and 
gauge invariance) can be parameterized using two functions F\ and F 2 . 



(f) = eu N (p') 



fi(Q 2 )/ + ^(Q 2 )^% 



u N (p) , (3) 



where Fi(Q 2 ) and F 2 (Q 2 ) are the Dirac and Pauli nucleon form factors respectively, k 
is the proton's anomalous magnetic moment, M is the proton mass, and q = p' — p, 
with Q 2 = —q 2 > 0. Employing this to evaluate the matrix element associated with the 
diagram in Fig. [TJ we find 



e 2 v 1 

M = ——u^p'^a, q]u e (p 2 )u N (p') 



\2\ 



un(p) , (4) 



where now q = p' 2 — Pi = p — p' is the four-momentum of the virtual photon that is 
exchanged. 

The spin- summed- and- averaged squared matrix element can be written as 



\ M \*= -rrL^H^—, (5) 



where the lepton and hadron tensors are both transverse with respect to the photon 
four-vector q, that is, 

q a H a(] = q^E afi = q a L aP = q^L a p = . (6) 

The lepton tensor can thus be replaced by 

L aP = Tr(p" 2 ^ Q i i> 2 fop) , (7) 

where terms proportional to the electron mass have been neglected, as they are suppressed 
by m e /m^. Straightforward evaluation then yields 

L a p = -4:q 2 (p' 2a P2f3 + P2 a P 2p ~ P 2 ■ V29 a fi) - 8p 2 ■ qp 2 ■ qg af3 . (8) 
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The evaluation of \M.\ 2 reveals that effects due to the Pauli form factor F 2 are suppressed 
by Q 2 /AM 2 . Below pion-production threshold this parameter is at most 0.02, and so 
in what follows we neglect the contribution to H a @ from the nucleon Lorentz structure 
ia^qy. For the proton, 

Fl (Q 2 ) = l-hr 2 p )Q 2 + 0(Q 4 ), (9) 

n 6 

with (r 2 ) 1 / 2 = 0.895(18) fm [5j, and so Fl (Q 2 ) = 1 up to a few per cent correction at the 
kinematics of interest here. Under these approximations, we obtain 

H aP = Tr((/ + M)7 Q (^+M) 7 /3 ) (10) 
= 4(p'V + p'V + (M 2 - p' ■ p)g aP ) . (11) 

Contraction of the tensors H and L then yields 

\M\ 2 = A 4 (g2)2 WW* ■ P')(P* ■ P) + 2 ^ ■ PW ■ P2) ~ 2M 2 (p' 2 ■ p 2 )} 

-4(p 2 • q)(p 2 ■ q)(p' ■ p) - 8MV 2 • ?)(P2 • q) + 8(p' • p)(p 2 • g)(p' 2 • ?)} ■ (12) 



Dropping terms which are suppressed by at least one power of Q 2 /M 2 relative to the 
dominant contribution, we are left with 

W? = ^Jch -M 2 -ml- Q%s - M 2 ) - M\m\ + Q>J + 0(Q ( 

where s = (p + p 2 ) 2 '■ 

Working now in the lab frame, and neglecting nucleon recoil, we have 



(13) 



E„ = E e - (14) 



^ = = m ^- 2E 2 + 2E e ^/E 2 - m£ cos0„ , (15) 

where 9^ is the angle between the outgoing muon and the incoming electron beam. The 
differential cross section is then 

,la Urr- f^T_ 1 ^ cos ^) ; (16) 



dCl^ A 4 



with E e = Ee/m^ and 



8E e A -6E 2 + 2E e JE 2 -l(l-4E 2 )cos9^ / x 

f(E e , cos 9,) = ^= -. (17) 



{2E e -2^E e -lE e cos6 ll -l) 



Integrating this over the muon solid angle we obtain 

E e 2 - 1 ( E e + \/E e 2 -l\ 

. (18) 



16na 2 v 2 H~2 
a = - — \ E P - 1 



A 4 v " 



E P 



In 
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For i? e just below pion threshold, the kinematic factor in the square brackets is about 2. 
Taking v ~ 200 GeV and A = 1 TeV results in a predicted cross section on the order of 
100 fb, which would definitely be observable. Including the nucleon-recoil terms neglected 
in the derivation of Eq. ffl8l) would result in corrections of order which are potentially 
as large as 20% or so, but do not change the order-of-magnitude of a. 

On the other hand, dimension-six BSM operators of the type in Eq. (J2j) do not induce 
effects mediated by low-momentum photons. If such operators do not contain additional 
derivatives they cannot produce powers of the nucleon mass in the numerator, and so the 
largest cross section they can yield is 

2 2 

am,, 

ff ~^> (19) 

which is suppressed by (^f-) 2 ^ 10~ 6 compared to the long-range mechanism depicted in 
Fig. [TJ Operators O that contain additional derivatives will be suppressed even further 
by at least one factor of the small parameter 

The result f fl8|) suggests that electron scattering from a proton target could provide 
access to BSM physics over a sizeable range of A. However, this prediction does not take 
into account the constraint on the /i — > coupling from the non-observance of this muon 
decay branch. As we shall see, this places stringent limits on the size of the cross section 
for the process in Fig. [TJ 

The operator in Eq. ([1]) produces an amplitude for the rare decay fi — > e-y that, in the 
muon rest frame, takes the form 

zev 

= -— ii e (-q) [ 7 °, j^K(0)^ , (20) 

where e is the photon polarization vector, and q is the photon four-momentum. From 
this, we obtain 

2e 2 v 2 m 3 , 



|AW 2 = ^r^- (2i) 

Converting this to a decay rate, and integrating over final electron states, we find 

2 S 

av 

r M ^e7 = ' (22) 



If we now use the result for the predominant muon decay mode 



G 2 F ml 

r M _ e ^p e = ig27r 3 , (23) 

(here Gp = 2 a/2 M t sin ig w * s the Fermi coupling constant, with Mw = 80.41 GeV the W- 
boson mass and 9w the Weinberg angle) we find that the branching ratio for /i — > e^f 
is 

BR(^e 7 )=9^4^. (24) 

But the factor ^ which appears here is the same as that in the pre-factor in Eq. (|T8|) . 
Consequently, we can eliminate this factor between Eqs. (fT8l) and (T241 to obtain 

a 1 

--^K^) 2 ^BR(^^e 7 ). (25) 

Using BR(/i e7) < 4.9 x 10 -11 Q, we see that 

a < 7.0 x KT 15 fb . (26) 

This is a model-independent constraint on the contribution to the cross section for ep — > /ip 
from photon exchange. Interpreted as a bound on A, we find A > 1.5 x 10 4 TeV. A similar 
bound on a was derived within the context of a MSSM calculation of the electron-nucleus 
to muon-nucleus cross section and BR(/i — > e 7 ) in Ref. [3J]. However, that number was 
seven orders of magnitude larger than the result of Eq. ( l26l) . Part of the difference arises 
from Ref. jsf's consideration of a target with Z = 70. 

One might ask why the apparent scale A in Eq. (CD) is so large — or, equivalently, why 
the coupling is "unnaturally" small. One possible explanation arises in the scenario known 
as "minimal-flavor violation" [sj]. There the physics beyond the Standard Model breaks 
the lepton-number symmetry of the Standard Model in the same fashion in which it is 
broken by neutrino mixing. This scenario can account for the small branching ratio for 
\i — > e 7 in a natural way as long as the product of the relevant neutrino mass and the 
scale of lepton-flavor violation A is smaller than v 2 . 

Regardless of what physics determines BR(/i — > ej), our calculations show that the 
bound on this quantity is sufficiently stringent to preclude the observation of any ep — ► \ip 
cross section from the diagram of Fig. [lj Indeed, the contribution of the operator in 
Eq. (CD) is constrained so strongly by the non-observation of this muon decay branch that 
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FIG. 2: Leading-order Feynman diagram for the process e~p — * [i~v^v e p. Solid lines represent 
leptons, the double line is a proton, and the dashed line is a W~ boson. The four-momentum 
carried by each external-state particle is indicated in parentheses. 

effects from diagrams with short-range operators of the form (j2J) are worth considering. 
In particular, if such EMU's involved a different scale A, with A < A of Eq. ([I]), they 
may produce a larger effect than fl26|) . However, the estimates provided above indicate 
that for A = 1 TeV, contributions from these dimension-six operators to the ep — > \ip 
cross section would be at most ~ 10~ 4 fb. The conclusion therefore is that beyond-the- 
Standard-Model physics is unlikely to result in any measurable production of muons when 
an electron beam impinges on a proton target. 



III. MUON PRODUCTION VIA STANDARD-MODEL ELECTROWEAK PRO- 
CESSES 



In this section, we evaluate the scattering cross section for the process e~p — > pTv^Vzp. 
The dominant contribution to the scattering amplitude comes from single-photon ex- 
change given by the Feynman diagram in Fig. [2j Applying the usual QED and Fermi- 
theory Feynman rules to the upper vertices in the diagram of Fig. [21 and using the 
single-nucleon current of Eq. (J3j) for the virtual-photon-nucleon vertex, we get 



iM 



ieu N (p')\<y"F 1 (Q 2 ) + K F 2 (Q 2 



2M 



G 



Ue(P2 ~ l)7a(l 



75j* , 2 
P2 



mi 



[V) 



For unpolarized electrons, the spin-summed-and-averaged squared matrix element can be 
expressed as 



\M\< 



2(Pa -m 2 e ) 2 ' 



= Tr [(/ + M) [^F x + ^V%F a } (p 1 + M) [fF x - ^^ e q e F 2 
L a/3 = Tr^ 27 a (l-75) ^/(l-Ts)] , 
W afll//3 = Tr [(/ 2 - Z)7 a (l - 75)(^2 + ^e)7^(^2 + m e )7„(/ 2 + ^e)7/?(l - 75)] (28) 



with I = lx + l 2 and 



(Q 2 



M 2 a 2 G 2 F 



Q 4 



The computation of the traces in W ailv p is facilitated by the Chisholm identity 



(29) 



(30) 



Performing the contractions, we obtain 

\M\ 2 = C j^2 7 U - I) ■ h] |4F 1 2 (A 1 ) - (^) 2 F 2 2 {A 2 ) - 2kF 1 F 2 (A 12 ) 
(A,) = 2(p' 2 ■ l 2 )[(p ■ p 2 )(p> ■ p' 2 ) + (p 1 ■ p 2 )(p ■ p' 2 ) - M\p 2 ■ p' 2 )} 

-p' 2 \p ■ p 2 )(p' ■ l 2 ) + (p' ■ p 2 )(p ■ l 2 ) - M 2 (p 2 ■ l 2 )\, 
(A 2 ) = 2(p' 2 -l 2 )[2(P-p 2 )(P-p , 2 )q 2 + 8M 2 (q-p 2 )(q-p' 2 ) + (p 2 -p' 2 )q 4 } 

-p' 2 [2(P ■ p 2 )(P ■ l 2 )q 2 + 8M 2 (q ■ p 2 )(q ■ l 2 ) + (p 2 ■ l 2 )q% 
(A 12 ) = 2(p' 2 ■ l 2 )[2(q ■ p 2 )(q ■ p' 2 ) + q 2 (p 2 ■ p' 2 )] - p' 2 [2{q ■ p 2 )(q ■ l 2 ) + q 2 (p 2 ■ l 2 )] , (31) 

where P = (p + p') and terms of 0(m 2 e ) have been dropped as p' 2 ^> m 2 . To eliminate 
the interference terms involving FiF 2 , we reexpress F\ and F 2 through the Sachs form 



factors 



to get 



G f 



Kq" 
AM 2 



F 2 and Gm = F± + kF 2 



(32) 



\M\ 2 
T 2 



1 CewM ~ • h\ J Ti _ o P2 ,7- 



P 2 ~ L P 2 • 

G T [P 2 (p 2 • / 2 ) - 2(P • p 2 )(P • l a )] + C 2 M [2(g • p 2 )(q ■ l 2 ) + (p 2 • Z 2 )g 2 ] 



G T [P 2 (p 2 • p' 2 ) - 2{P ■ p 2 ){P ■ p' 2 )] + G 2 M [2(q ■ p 2 )(q ■ p' 2 ) + (p 2 ■ p' 2 )q 2 ] ,(33) 
where G T = {G% + tG 2 m )/{1 + r) with r = Q 2 /AM 2 . 
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In the laboratory frame, the differential cross section is given by 
1 f d 3 k d% d 3 p' d 3 l 



(a ' AMEj (2nf2E^ (2^2E^~- (2nf2E p , (2nf2E Ue 1 1 

x(2tt) 4 / d 4 p , 2 5^(p 2 +p-p' 2 -p') [ d*l5W(l-h-l 2 )5M(j/ 2 -l-l Ve ),(M) 



where E e is the energy of the incoming electron in the laboratory frame, l Ve is the electron- 
neutrino's 4-momentum and the last two delta functions and integrals have been inserted 
as unity to simplify calculations. The steps given in Appendix El then allow us to obtain 
from Eq. an expression for the differential cross section per unit solid angle subtended 
by the detected muon at fixed beam energy: 



dQ, 32^M^E e J m ^ " »J qh ql 



dE^E^-ml -fdql + 2Mq d(cos6 q ) 



q J cos 



30' 



x / d(j) q (p 2 -q- h, 



/4"^2 — ——To %. 

P 2 2p' 2 



(35) 

(P2=(P2-</),g 2 =-2Af<3o) 



J 2 

The limits on the go integral are determined by the electron beam energy. The lower 
limit of the integral arises because energy transfer to the proton without any momentum 
transfer is not possible in elastic scattering: the target recoils. In the limit that the muon 
is produced at rest (I2 = 0), an analytic expression for both the upper and lower limit 
of the go integration can be obtained. This can guide intuition on the importance of 
collective excitations when the target is replaced by a heavy nucleus (see Sec. |V|). We find 



± _BtVB 2 -AAC 
% ~ 2A 5 

A = 4(E e + M - m M ) 2 - 4E 2 e ; C = m\{2E e - m 



B = 2 [2{E e + M- mJm^Ee - m„) - AE 2 e M] . (36) 



Furthermore, 



2 

Tn 

% = —, r <E e -E u , (37) 

y ° 2(M-m M ) M V ; 

so that the minimum electron beam energy E™ n for muon production is then determined 
by requiring q$ = E^ m — m M , which yields 

Ef* = - m/j) =111.6 MeV. (38) 

The quantity q$ is thus at least 6 MeV, whereas 

+ _ (2E e - m M ) 2 
% ~ 2(2E e — m M + M) (39) 
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which is always less than E e — m M . The corresponding Q 2 ranges from 0.01 GeV 2 to a 
maximum of 2Mq$ . The angle between the electron and muon neutrinos (whose masses 
are neglected) is constrained to < Z(li,l Ve ) < tt/2 by the step-function Q(Q 2 ). The 
maximum energy in neutrinos is £'™ ax = E e — q$ — m M = E e — E™ m . 

It is noteworthy that the differential cross section is independent of the azimuthal 
angle M . Dependence on M appears explicitly in the matrix element through the dot 
products (P ■ Z 2 ) and (q ■ l 2 ) as well as implicitly in the step function 9 {{k — q — h) 2 ), 
but this dependence drops out once the d<p q integration is performed. Only differences of 
azimuthal angles {<p q — 0^) appear in the d<p q integrand as well as in the limits on this 
integral, so the integral remains invariant. Therefore experiments to measure this process 
are characterized by M alone, and the 0-independence can be used to increase the total 
number of counts (thereby decreasing the statistical error) by positioning several detectors 
in an annulus at the same 9^. 

The integrals in Eq. ( l35l) were evaluated numerically, details of which are presented in 
Appendix [HI The result for the differential cross section da/dVl^ as a function of electron 
beam energy at a fixed value of # M = 7r/3 will be presented in Sec. IIVI The corresponding 
total muon-pro duct ion cross section rises from 3 x 10~ 16 fb at E e = 120 MeV to 4 x 10~ 13 
fb at E e = 140 MeV. 

The cross section for the reaction e~p — > [i~pi> e i?^ is thus much smaller than the low- 
energy approximation to the Rosenbluth cross section for ep elastic scattering flo| : 



At E e = 140 MeV this is ~ 10 7 fb for all but forward angles where the Coulomb singularity 
occurs. Thus, muon production via standard-model electroweak processes is down by 20 
orders of magnitude as compared to elastic electron-proton scattering. Much of this 
suppression comes from the extra factor G 2 F l 2 (1^ ■ l v ) (with l x the four-momentum of 
lepton x). Numerically, G% ~ lO" 10 GeV" 4 , l Ue ~ 0.01 GeV, Z M ~ 0.1 GeV, and l Vfi ~ 0.01 
GeV, and so this factor already implies a suppression of 17 orders of magnitude. 

Further suppression occurs due to the lower bound on the virtuality of the exchanged 
photon, which was explained above, and persists even if the muon is detected at 9^ = 0. 
The Coulomb divergence that is manifest in Eq. (HUj) as 9 — > and is associated with 
Q 2 — > does not appear in e~p — > ^pu e u^. The suppression relative to the Rosenbluth 





(40) 
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cross section is therefore even more severe than is implied by the dimensional analysis in 
the previous paragraph. The lack of enhancement from the exchanged photon going soft 
is an important feature of muon production. 

The cross section predicted for standard-model electroweak [T~ production at the 
largest energy considered here, E e »s m^, is thus two orders of magnitude larger than 
the largest BSM cross section predicted by Eq. (126]) . We limit the indicent electron en- 
ergy to E e < m n as for incident energies larger than the pion mass, strong-interaction 
processes involving the production of an on-shell pion which then decays to a muon will 
swamp the purely electroweak diagram of Fig. [2j In the following section, we show that 
a sub-threshold version of the pion-production process yields muon-production cross sec- 
tions that are significantly larger than those obtained through the mechanism discussed 
in this section. 



IV. MUON PRODUCTION VIA SUB-THRESHOLD PION PRODUCTION 



In electron-proton collisions, muons can also be generated through processes in which 
a virtual photon couples to a pion which then decays into a muon and a neutrino (Fig. 
[3]). This is a manifestation of the "pion cloud" of the nucleon, and the processes depicted 
in Fig. [3] are possible even when the pion in the intermediate state is virtual, i.e. E e is 
significantly below the pion mass. Whereas only negatively charged muons can be pro- 
duced in the processes considered thus far, this strong-interaction process yields positively 
charged muons, together with a neutron in the final state. In this section, we evaluate 
the differential cross section for the reaction: ep — > enji + u tM . 

Our calculations are performed using chiral perturbation theory (xPT), the low-energy 
effective field theory of QCD. incorporates QCD's (broken) SU{2)l x SU{2)r sym- 
metryas well as the pattern of chiral-symmetry breaking in QCD (for a recent review see 



Ref. 



111]). Reactions involving pions, nucleons, and photons (either real or virtual) can be 



straightforwardly and systematically evaluated using %PT, as long as the energies involved 
are well below the excitation energy of the A(1232). Here, we perform a tree-level calcula- 
tion of the process of interest using the leading-order %PT Lagrangian. For diagrams (1) 
to (3) in Fig. [3] our calculation is equivalent to evaluating the amplitude for charged-pion 
electroproduction at leading order 0(e). Such a leading-order calculation is known to 
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FIG. 3: Leading-order xPT contributions to ^-production in electron-proton scattering. The 
double line denotes the nucleon, the thin solid lines denote leptons and the dashed lines denote 
pions. Particle momenta are indicated in parentheses. The numbers in parentheses correspond 
to the individual amplitudes computed below. 

give a reasonable description of the available data for charged-pion photoproduction near 
threshold [l2] . 

The leading-order xPT Lagrangian describing the interactions between pions, photons, 



and nucleons is given by 



(2) 

Here CXI denotes the leading-order Goldstone-boson Lagrangian 

£™ = S-Tr[D^U(D»Uy} + &Tr( X ^ + Ux ] ) , 
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C 



Nir ~ *—-n-7T • 



X PT - J- NlT 



(41) 



(42) 



where to leading order in quark masses the matrix x is m\ times the identity matrix, and 
denotes the lowest-order Lagrangian involving baryons: 



9a 



(43) 



In the above equations, the pion fields are collected in the matrix U = exp(zr ■ 7c/f n ), 



whereas the fields u = exp(ir -tt/ (2f n )) and = i( u 



o 

'd^u — ud^w). The quantities 9 a and 
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M denote the axial coupling constant and the nucleon mass, respectively, in the chiral 
limit. The covariant derivative acting on the pion matrix is defined as 

D^U = d„U - ir„U + iUlp , (44) 

where r M and Z M denote the appropriate external fields and external electromagnetic fields 
Afj, are coupled to the pion field by setting r M = Z M = A^. The (chiral and U(l) em ) 
covariant derivative acting on the nucleon field is then 



(45) 



As our leading-order computation involves only tree-level diagrams, we can employ the 
relativistic Lagrangian for nucleon and pion fields without concern about contributions 
from loop graphs that might violate power counting [ijj . We will therefore use relativistic 
Feynman rules in what follows. The amplitude corresponding to each of the diagrams 
contributing to /i + -production at leading order (see Fig. [3]) is then 

M.\ = —iV^VudGpf-^g^NNe 2 — — ■ — r- -— -— 

[K + q) z — m z k z — m z q z 

Xu v {k + q - p 3 )(# + 4){l - 7 5 )^(Ps) 
XM e (p 2 - q)(2# + <i)u e (p 2 ) 

xu N (p - k)^ 5 u N (p) , (46) 
M 2 = -iV^VudGpUgirNN e° 



xu v (k + q - p 3 )(# + 4)(l - 7 5 )^(P3) 



(k + q) 2 — m 2 (p + q) 2 — M 2 q 2 

'm( 

xu e (p 2 - q)7 M M e(P2) 

xujv(p - k) 7 5 (^+ 4 + M) 7 ^u 7V (p) , (47) 



M 3 = -iV2V ud G F f n g nNN e 2 -^- 



1 



2M (k + q) 2 - q 2 
xu v (k + q-p 3 )W + 4){l - 7 5 K(Ps) 

XU e (P2 - q) 7 M «e(P2) 

XM 7V (p-k)7 5 7 M M 7V (p) , (4* 
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FIG. 4: The differential cross section in femtobarns per steradian versus the energy of the 
incident electron beam. The solid line gives the result for production of a /i~ via the process 
e P ~^ ^~P^nV e (see Sec. IIIip . and the dashed line is the result for production of a /u + via the 
reaction e~p — ► e~n/i + i^ (see Sec. lIVp . Both results were evaluated for a representative muon 
angle # M = 7r/3. 



_M 4 = -iV2V ud G F f^g v 



1 



nn e 



G°3 ~~ ?) 2 ~~ m u ^ 2 ~~ 777,2 ? 2 



x«,(k + q - p 3 )(# + tf)(l - 7 5 )(]/3- + «v)7m<V(P3) 
xw e (p 2 - q) 7 M M e (p 2 ) 

XU N (p - k)7 5 MAr(p) , 



(49) 



where g nNN = ™ at this order, and we adopt # A = 1.26, / w = 92 MeV, M = 939 MeV. 
Note that the crossed counterpart of diagram (2) is zero if only leading-order couplings 
are considered, as this process involves a neutron in the final state. 

We evaluate the matrix elements M.\— M.a using the package FeynCalc [ijj]. This 



produces an expression for the spin-averaged-and-summed squared matrix element \M.\ 2 
that is lengthy and not particularly illuminating. The differential cross section is then 

(27r) 4 «5(p' + p' 2 + p 3 + p v - p 2 - p) \M | 2 , (50) 



rffi^ 2M 2£ 2 J (2tt) 12 2E; 2£ 2 2E 3 2E v 



16 



where p u = (E v , p u ) is the four-momentum of the outgoing neutrino, and the four vectors 
p' = p — k, p' 2 = P2 — Q, and p% (which is the outgoing muon momentum) are written 
in a similar fashion. We evaluate the integrals in Eq. (1501) by Monte Carlo integration, 
obtaining a result that is numerically stable to better than 5% accuracy. 

The results of our calculation are shown by the dashed line in Fig. HJ where the energy 
dependence of the differential cross section at a representative outgoing muon angle of 
Op = 7r/3 is displayed. The solid curve in this figure shows the differential cross section for 
the production of pT through the photon and W~ mediated mechanism of the previous 
section. The dashed curve shows results for the production of /i + 's through virtual-photon 
exchange discussed in this section. The differential cross section for p + production is four 
to five orders of magnitudes larger than that for pT production. Even the /i + -production 
cross section is, however, very small: of order 10~ 9 fb at the largest energy considered 
(E e = 140 MeV). The variation of the cross section with the angle 9^ is one order of 
magnitude for both cross sections, so we predict a total cross section for p + production 
of order 10~ 8 fb just below the pion threshold. 

The dependence on energy of the total cross section for the processes e~p — > e~np + v ll 
and e~p — > pTpv^Vf, (see Sec. IHIj) is shown in Fig. |5j Also shown in Fig. [5] is the bound of 
Eq. (126]) for ep — > pp from photon exchange. Note that in Fig. [5] we do not display results 
for energies exceeding 140 MeV, because above that energy the nir + channel opens and 
p + 's are copiously produced through the decay of real pions. 

Even at E e = 140 MeV, the cross section for p + production via strong interactions 
is many orders of magnitude larger than the cross section for BSM muon conversion in 
Eq. (|26|) . Indeed, it may be competitive with BSM mechanisms even if dimension-six BSM 
operators that induce EMU are not suppressed by, e.g. minimal lepton-flavor violation. 
Thus, any experiment that searches for EMU on a proton target via BSM processes 
should discriminate between the desired reaction and the channel ep — > enp + u fJi . Such a 
discrimination requires either detecting the outgoing electron, or detecting the charge of 
the final-state muon. 
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FIG. 5: Total cross sections in femtobarns plotted against the energy of the incident electron 
beam. The solid line gives the result for production of a /i~ via the process e~p — * /i - pi^z^ e (see 
Sec. IIIip . and the dashed line is the result for production of a fi + via the reaction e~p — > e~np + v^ 
(see Sec. HVj) . The dotted line is the bound on BSM contributions obtained for ep — > fip in 
Section [II] by considering the dimension- five BSM operator and the non-observation of the decay 
H ery. 

V. MUON PRODUCTION IN ELECTRON-NUCLEUS SCATTERING 

If the proton target were replaced by a neutron target (e.g. via the use of neu- 
trons bound inside a deuterium nucleus), the incoming electron can also interact with 
the neutron through its magnetic moment. From Eq. (|27I) . the matrix element for 
magnetic-moment interactions introduces an extra factor — relative to the dominant charge 
interaction — of q^/M. This translates to a factor Q 2 /M 2 « q^/M m 0.01 in the cross sec- 
tion. Thus, the electroweak process en — > /i~z/ e z/ M n yields a smaller muon-production 
cross section than in the case of a proton. In contrast, the process ep —>■ enfi + u^ dis- 
cussed in Sec. HVl is associated with an isovector matrix element at leading order in xPT, 
and so the cross section for production of muons will be as large for a neutron target as 
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for a proton target. In this case the reaction is, however, en — > epfi^u^. For neutrons 
this muon-pro duct ion reaction provides a signal that cannot be distinguished from BSM 
electron-muon conversion by detection of the charge of the final-state muon. (As an aside, 
we note that in the case of Z° exchange, scattering off a neutron is more favorable due 
to its much larger weak charge as compared to the proton. However, the appearance of 
an extra factor of Gf renders the cross section due to electromagnetic- weak interference 
terms negligible when compared to photon exchange.) 

The cross section for muon production is enhanced when the electron scatters off a 
heavy nucleus. As an illustrative example, and to estimate the expected enhancement 
over the nucleonic case, we consider electron scattering on a lead nucleus ( 208 Pb). The 
low energy and low three-momentum transfer region is probed in conventional nuclear 
spectroscopy. In this region, the elastic peak appears first, although at q^ h = — g 2 /2mpb 
instead of q^ 10 ' = —q 2 /2M, so that the exchanged photon appears to be two orders of 
magnitude softer. In fact, though, the requirement of producing the muon implies that 
the photon's virtuality is unchanged from Q 2 ~ m 2 ; as the differential cross section scales 
as ~ 1/(Q 2 ) 2 = l/(2mpbgo b ) 2 considering only the effect of the heavier target yields no 
particular enhancement. This fact can also be verified by counting powers of the target 
mass and energy transfer in the expression for the differential cross section in Eq. (135]) . 

The large charge of lead (Z = 82) does tend to increase the cross section, although 
nuclear elastic form factors offset this effect substantially. The typical three-momentum 
transfer involved in elastic scattering is |q| = \/q 2 + 2mpt,go ~ m ti, which corresponds to 
a spatial resolution of about 2 fm. But Rpt, ~ 7 fm is the typical size of the charge dis- 
tribution in 208 Pb as determined by fitting a conventional 2-parameter Fermi distribution 



for a spherical nucleus .16], so we do not expect the lead nucleus to respond coherently 
to the electromagnetic probe. This can be quantified if we approximate the elastic form 
factor by the diffraction pattern from a spherical charge distribution of radius -Rpb- In so 
doing we obtain an overall factor relative to the proton case of: 



Z 2 F 2 (\q\) ss Z 2 



3ji(| 


q 


Rpb) 




q 


-Rpb 



(51) 



where ji is the spherical Bessel function of the first kind. For lead, F 2 ~ 10 2 at the 
Q 2, s of interest here, in good agreement with form factors extracted from data on elastic 



scattering from 208 Pb in this region of energy and momentum transfer 
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17| . Therefore, if 



we replace the target proton by a lead nucleus, we expect an overall increase of the cross 
section by Z 2 F 2 ~ 82 2 x KT 2 « 67. 

Elastic scattering is not the complete story, however, because the maximum energy of 
the exchanged photon is ~ (E e — m /t ) ~ 30 MeV, which is sufficient to excite a tower 
of collective states. Studies of inelastic form factors of the first few excited states (for 



example, the 3 octupo 



to the elastic peak 
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e in Pb) reveal a suppression of ~ 10 or more compared 



19]. The width of these excited states is also small (0.1 MeV 
for the 3~ state), therefore, at low energies of the exchanged photon, the contribution of 
the elastic peak is dominant. At slightly higher energies, go ~ 10 MeV, giant monopole 
and multipole resonances can be excited. These resonances are of empirical importance 
in studies of nuclei as they carry non-zero isospin. Although the resonances have large 
widths (1-5 MeV), their contribution to the cross section will also be smaller than that 
from the elastic peak. 

Finally, we inquire whether quasi-elastic scattering should be taken into account. By 
quasi-elastic scattering, we are referring to those events in which a muon is produced 
and a nucleon is knocked out of the nucleus. This phenomenon requires the additional 
kinematic restriction that the three-momentum transfer exceed the Fermi momentum of 
the nucleon in the nucleus. Therefore, the relevant energy regime is now defined by the 
conditions 0(|q| — kp) and the theta functions imposed above, i.e, Q(E e — g — E^) and 
©((& — q — h) 2 )- These two theta functions are unchanged from the nucleonic case as they 
originate from the kinematics of the leptonic portion of the process, which is unaffected 
by changing the target from a nucleon to a nucleus. These restrictions imply that the 
maximum value of (k — q — Z 2 ) 2 is given by 

(E e - g ) 2 - E 2 e -k 2 F + ml~ 2m^E e - g ) + 2E e k F cos(kq) > , (52) 

where the inequality imposed by the theta function 6((/c — g — Z 2 ) 2 ) is satisfied so long 
as < g < g^, where g^ is the lesser root of the above quadratic in g . Clearly, this 
requires that g^ > 0, which is equivalent to the condition 



m. k% — m 2 
C °^^t + ^Ejf- (53) 



As |cos(A;g)| < 1, we obtain the restriction 



Ee > |'*£±^'] . (54) 
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If we assume a simple picture of the nucleus with constant density p ~ p nuc = 0.16 
fm -3 , then kp « 260 MeV, which implies that E e > 187 MeV. This exceeds the pion- 
production threshold in ordinary electron-nucleus scattering (no muon production). It 
is highly desirable that the electron beam energy not be above the pion threshold at 
around 140 MeV, and in this case we need not include the contribution from quasi-elastic 
scattering, since Eq. (|54|) makes clear that it is important only at energies well above pion 
threshold. This is significantly different to the usual situation in inelastic electron-nucleus 
scattering (i.e. without muon production), in which pion production occurs at energies 
that exceed the quasielastic peak. When muon production happens, additional kinematic 
restrictions (viz., the energy cost of producing a muon) imply that the quasielastic peak is 
only important at energies that exceed the threshold for pion production. This is another 
distinguishing feature of the muon-pro duct ion process. 



VI. CONCLUSIONS 



We have examined the possibility of discovering physics beyond the Standard Model 
through lepton-flavor violation in fixed-target electron scattering. Our main findings can 
be summarized as: 

• We have obtained a model-independent constraint on the magnitude of LFV in 
electron-nucleon scattering from beyond-the-Standard-Model effects using a general 
low-energy effective interaction with couplings constrained by experimental bounds 
on the nonobservance of // — > ej. The cross section for LFV from the lowest- 
dimension operator, a < 7 x 10~ 15 fb, is too small to be experimentally accessible 
with current technologies. The contribution of higher-dimension LFV operators to 
ep — > \ip could be larger, but is still unobservably small at present. This is in accord 
with similar estimates that have been made previously within specific extensions of 
the Standard Model UM- 

• We have identified two main sources of background in inclusive ep scattering within 
the Standard Model when only the energy of the outgoing muon is measured, and 
performed detailed calculations of the relevant cross sections. The reaction e~p — > 
jjTVePyj) is the principal background if the charge of the muon is measured, and its 
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cross section varies from the order of 10 16 fb at incident electron energy E e = 120 
MeV to 10- 13 fb at E e = 140 MeV. 

• If the charge of the muon is not measured, the dominant source of background comes 
from /i + s produced by the decay of virtual pions. Leading-order chiral perturbation 
theory gives this reaction's total cross section a(ep — > en/i + z/ M ) to be about 10~ n 
fb at incident electron energy E e = 120 MeV and ~ 10 -8 fb near the pion thresh- 
old. This background swamps any LFV signal in ep scattering unless the outgoing 
electron is also detected or/and /i + events are vetoed. 

• Using a heavy nucleus as a target enhances both the desired LFV effects and the 
background. At the low energies carried by the exchanged photon in e~p — > {i~pv e v^ 
the role of collective nuclear excitations can be neglected in comparison to the 
leading effects of elastic scattering from a finite-size target. This could enhance the 
cross section for e~p — > fi^puei?^ by as much as two orders of magnitude, but the 
cross section is still too small to be experimentally detectable at present. 



APPENDIX A: PHASE-SPACE EVALUATION FOR FINAL STATE p^v e 

In this appendix, we explain how to obtain Eq. fl35l) from Eq. fl34|) . We first employ a 
useful relation for elastic scattering 

f d 3 p' 1 / q 2 \ 

J 2e/ iP + q ~ P>) = 2M 6 \2M + % ) Wh6re % = P ™ ~ P2 ° ' (A1) 
which enables Eq. (1341) to be rewritten as 

' ' "■ 8^^-h-l Ve -l 2 )\M\\ (A2) 



J (2^)32^(2^)32^ 

For massless neutrinos, the phase-space integrals over neutrino momenta in the second 
line of Eq. (|A2[) can be rewritten as 

/ <^ W <i-Pi + « / (2 ^k(2^/ ' 4,(L -''-'- )W - <A3) 
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Noting from Eq. (I3"3l) that Eq. (1A3|) has a factor (p' 2 — l)-l\ = l Ve ■ h, the integrals over 
d 3 l Ue and d 3 h are |20j 



d 3 Zi 



d 3 Z, 



-5^{L-h-l Ve ) l Ue .h 



7TL 2 



(2tt) 3 2^ (27r)32^ e v 1 ^ 4(2tt 

Using the resultant expression for Eq. (IA3I) in Eq 



fl(L )fl(L 2 ) , (A4) 
2]), performing 



where L 2 

the <i 4 L and d 4 p' 2 integrations with the aid of corresponding delta functions, and using 
Eq. (|29|) . we obtain 



da 



a 2 G 2 F 



Att 2 M 2 E p / 2E. 



d% 



d*q 



5 



T 



2M 



+ Qo 



)l(q,k 



(p 2 -q- h) 



P 2 -h 

I 4 
. P 2 



e ((p 2 - q - l 2 f) Q(E - Ep- - g ) (A5) 

p' 2 =(p2-q) 



where T±,T 2 are given by Eq. fl33p .With the aid of the only remaining delta function, the 
J d 4 q can be recast as 



d 4 q 



M 



2M 
dq 



+ qo i{q, k) 



q 2 + 2Mq Q(E e -Ep- - q ) 



x 



-2q M) 2 

J dn q e ((p 2 - q - l 2 ) 2 ) (p 2 -q- h) 



.P' 2 



V 2 



(A6) 



(p' 2 = (P2-q),q 2 =-2Mq ) 



The step functions Q(E — q — E^ 



and ((p2 — q — h) 2 ) provide the upper and lower 



limits on the dqo integral. The latter step-function also provides bounds on the angular 
integrations involving dcos9 q , d<p q . This determines the support for the various integrals 



as [?o'9o]> [cos0' cos0"], 



and leads to (1331) . In our numerical calculations, we have 



used a constant value for Gp(Q 2 = 0.01GeV 2 )=1.05 x 10 5 GeV 2 as determined by its 



Standard-Model running in the MS scheme [21 ]. 



APPENDIX B: NUMERICAL NOTES 



The integrals in Eq. (1351) are performed as follows. We choose the +z axis to be along 
the electron beam direction. The polar angle 9^ is measured from the +5-axis in the 
vertical plane containing this axis. The azimuthal angle M is measured anti-clockwise 
from the (arbitrary) — z axis in a plane containing this axis. The position of the detected 
muon is then uniquely specified by the angles 6^ and Once the position of the muon 
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is specified as above, for a fixed momentum we can determine the range of go for which 
the step functions in Eq. fl35|) do not vanish. This procedure determines the bounds on 
cos9 q at fixed beam energy E e , from which bounds on (p q follow. The numerical evaluation 
of the multiple integral is then performed using standard quadrature methods. At the 
low Q 2 values involved here the Q 2 -dependence of Ge and Gm induces a correction of 
5-10% in the cross section, as compared to using their q 2 = values. We have taken 
this into account in the numerical results presented in Sec. Hilt using a standard dipole 
parameterization obtained from studies of e~p scattering: 

G E (Q 2 ) = ^^ = j . (Bl) 

1 + k (l + Q 2 /0.71GeV 2 ) 2 1 ; 

For the range of Q 2 relevant to the process considered here, this parameterized form 
is accurate to better than 1%. 
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